Chandra Contaminant Migration Model 

Doug Swartz & Steve O Dell 

Motivation: 

The accelerated accumulation observed since -2012 cannot be reproduced 
under the old assumptions of one (or more) gradually depleting 
contaminant source(s) — any new paradigm must consider source(s) with 
rates increasing with time (but consistent within S/C trends). 

Methods: (the usual) 

(1) Simulate migration (vaporization/deposition) of potential 
contaminant(s) by solving a set of 1 st order ODEs through explicit 
(forward) difference scheme and within thermal & geometric model 
constraints. 

(2) Compare time evolution of material mass column density to 
observations at select locations on OBF filters to identify best candidate 
contaminant properties. 
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Chandra Contaminant Migration Model 

Simulation Constraints: 

♦ Spacecraft geometry (surface areas & view factors) and time-averaged 
surface temperatures as determined by Neil Tice (LM) using flight data 
inputs to Thermal Desktop™ radcad (geometric) & sinda fluint 
(radiative/conductive thermal) modules; accounts for effects of DH heaters 
being turned off -2008.4 (t~9 yr) 

♦ X-ray observations of contamination effects throughout the mission as 
documented in CALDB contamination model N0009 (June 2014); 
constrains time evolution and spatial distribution of contaminant optical 

depth, r(t,x,y), proportional to the mass surface density, (l (t,x,y) 

♦ Temperature dependence of vaporization rate follows Clausius- 

Clapeyron relation p(T) = p(To) (T/To) 1/2 exp(-AH(l/T-l/To)/R) 

where p(T) is the mass vaporization rate at temperature T, AH is the 
vaporization enthalpy, and R the universal gas constant 
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Chandra Contaminant Migration Model 

Simulation Parameters: 

♦ Material (enthalpy, AH, and reference temperature, To , fixed); 
vaporization rate normalization p(T 9 ) is a free parameter 

♦ Declining source rate : S(t)=A(l+B*exp(-t/C)) where A, B & C are 
chosen to follow the early build-up (t<2008) at the center of S3 

♦ Rising source rate : S'(t) a A'exp(-B'(l/T-l/To) where the normalization, 
A' is a free parameter and B' is given by the enthalpy AH of the supposed 
contaminant. Physically, A' includes a multiplicative scale factor 
representing the outgassing ‘escape efficiency’. 
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Chandra Contaminant Migration Model 


ToBA(t) = To [1+0.006 exp((t-8.1)/4.0) ] 


OBA CONE AVERAGE TEMPERATURE (OBACAVE) [50 0 ♦/- 2 5], TELESCOPE, (DEG F) 
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...the Clausius-Clapeyron equation for a time-dependent TobA 
( AH and To are those of the contaminating material) 
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Chandra Contaminant Migration Model 
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Chandra Contaminant Migration Model 

SUMMARY 

♦ Volatility of contaminants can be tightly constrained by 
simulations. 

♦ Accelerated buildup since 2010 has 30* higher volatility at 
-60°C ; likely to ‘clean’ readily at elevated temperatures. 

♦ Current exponentially increasing trend in OBA temperature is 
predicted to lead to extremely rapid buildup of ‘second’ 

contaminant (5* source rate for 10°C increase in 3 yrs; 30* in 5 yr; 

but polynomial fit implies only 3* in 5 yr). 

♦ At t~16 yrs, 0.4 grams of contaminant is within the S/C with 

only 5% (20mg) on OBFs (-270 //g/cm 2 ). 
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Optical Bench 
Closeout 


Geometry Model 


Collimator 



Camera Body 

(ACIS Housing) 


From NeilTice/LMC 

Thermal Desktop (finite element): 
RADCAD 

to calculate geometric view factors 

SINDA FLUINT 
to calculate temperatures 
738 nodes, 121 OBF-1,203 OBF-S 

1 86x 1 86p/node (I), 1 46x2 1 2p/node (S) 


Snoot 
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Chandra Contaminant Migration Model 

Results 

Mass Column of Octadecane (C 18 H 38 ) at t=9 years 



low volatility (0. 10) 

“deposition” dominated: central 
regions have highest accumulation 
because center views more nearby 
cold surfaces, pattern is asymmetric 


high volatility (2.50) 

“thermal” (vaporization) dominated: warm 
central regions begin to clean, pattern 
follows local temperature distribution with 
more material near cold edges 
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